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Abstract 

The neutrino emissivities resulting from direct URCA processes in neutron 
stars are calculated in a relativistic Dirac-Hartree approach in presence of a 
magnetic field. In a quark or a hyperon matter environment, the emissivity 
due to nucleon direct URCA processes is suppressed relative to that from 
pure nuclear matter. In all the cases studied, the magnetic field enhances 
emissivity compared to the field-free cases. 

PACS numbers: 26.60.+C, 21.65.+f, 12.39.Ba, 97.60.Jd 
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Typeset using REVT^ 



Neutron stars are born in the aftermath of supernova explosions with interior tempera- 
tures T ~ 10^^ K, but cool rapidly in a few seconds by predominant neutrino emission |l| 
to T < 10^° K. Neutrino cooling then dominates and lasts for t ~ 10^ — 10^ yr and sub- 
sequently photon emission takes over when T ~ 10*^ K. Since the long term cooling of the 
young neutron stars (T ~ 10^ — 10^° K) proceeds via emission of neutrinos primarily from 
matter at supranuclear densities within the core, the study of the cooling of neutron stars 
by examination of neutrino emissivities may provide considerable insight into their interior 
structure and composition. 

For a long time the dominant neutrino cooling mechanism has been the so-called standard 
model based on the modified URCA processes 

{n, p) + n ^ {n,p) + p + + Fg, 
{n,p) + p + e~ {n,p) + n + Ue. (1) 

The ROSAT detection of thermal emission from neutron stars indicates the necessity of 
faster cooling mechanism in some young neutron stars, in particular the Vela pulsar. Faster 
neutrino emission than the standard model was proposed by invoking pion or kaon 
condensates which have neutrino emissivities comparable to that from the /3-decay of quarks 
in quark matter (consisting of u, d and s quarks) 

d ^ u + e~ +17^, u + e~ ^ d + Ue- (2) 

A similar relation for s-quark /3-decay may occur and is obtained by replacing d by s quark 
in Eq. (2). The most powerful energy losses, expected to date, are produced by the so-called 
direct URCA mechanism involving nucleons |^ 

n —>■ p + e~ +T'e, p + e~ ^ n + Ug. (3) 

The threshold density of this process is however considerably larger than that of the modified 
URCA process. 

Observations of pulsars predict large surface magnetic field of Bm ~ 10^^ G @. In the 
core the field may be considerably amplified due to flux conservation from the original weak 
field of the progenitor during its core collapse. In fact, the scalar virial theorem predicts 



large interior field ~ 10^^ G or more |]rT|, and these fields are frozen in the highly 



conducting core. It has been demonstrated [11 1 that when the field is comparable to 



or above a critical field -B^\ the energy of a charged particle changes significantly in the 
quantum limit; the quantum effects are most pronounced when the particle moves in the 
lowest Landau level. The phase space modifications stemming from the strong magnetic 
field in the core are expected to influence the neutrino emission rate from young neutron 
stars. 

In this communication we evaluate the neutrino emissivity for the nucleon direct URCA 
process of Eq. (3) in presence of a magnetic field Bm, and demonstrate that it would lead to 
more rapid cooling in the core. (A straightforward extension of the emissivity for nucleons 
into the quark sector may also be obtained in a magnetic field.) For this purpose, we 
consider a npe matter in /3-equilibrium within a relativistic Dirac-Hartree approach in the 



linear a-co-p model ||T2 . 
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At the neutron star core at temperatures well below the typical Fermi temperature of 
Tp ~ 10^^ K, the nucleons and electrons participating in neutrino producing processes are all 
degenerate (the and z/g are free) and have their momenta close to the Fermi momenta pp^, 
where i = n,p, e. Since neutrino and antineutrino momenta are ~ kT/c -C pPi, the nucleon 
direct URCA process is allowed by the momentum conservation when PFp+PF^ ^ Pf,,- Since 
matter is very close to /5-equilibrium, the chemical potentials of the constituents satisfy the 
condition /i„ = fip + fie- (Henceforth we set ^ = c = = 1.) 

Employing the Weinberg-Salam theory for weak interactions, the interaction Lagrangian 
density for the charged current reaction (3) may be expressed as C^^^ = {Gp/ V^) cosOJ^j^^, 
where Gp ~ 1.435 x 10^^^ erg cm~'^ is the Fermi weak coupling constant and the Cabibbo 
angle. The lepton and nucleon charged weak currents are respectively, = ■?/'47^(l — J5)ip2 
and = ip^'^'^{gv — gAl5)4'i- Here, and in other formulae to follow, the indices i = 1 — 4 
refer to the n, Ve, p and e, respectively. The vector and axial-vector coupling constants are 
gv = 1 and gA = 1.226. 

The emissivity due to the antineutrino emission process in presence of a uniform magnetic 
field Bjn along z-axis when both the electrons and protons are Landau quantized is given by 



f Vd^pi 


f V<fp2 


1 (27r)3 J 


' (27r)3 J- 



X 



qB„^U/2 2tT J-pFp 27r 2tT J -pPe 27r 

??max ^max 

E Y: E^WfJMil - /(p3)][l - /(P4)] , (4) 



where ?7max and rj'j^i^^ are respectively the maximum number of Landau levels populated for 
protons and electrons. The prefactor 2 takes into account the neutron spin degeneracy. The 
Pi = {Ei, Pi) are the 4- momenta and E2 the antineutrino energy. The functions f{Ei) denote 
the Fermi-Dirac functions for the ith particle. The transition rate per unit volume due to 
the antineutrino emission process may be derived from Fermi's golden rule and is given by 
Wfi = {\J^ fil"^) / (tV) . Here t represents time and V = L^LyL^ the normalization volume. 

is the squared matrix element and the symbol (■) denotes an averaging over initial 
spins and a sum over final spins. The matrix element for the V — A interaction is given by 



Mfi = -^J d'x {gv - gAi,) MX) (i - 75) MX) • (5) 

In presence of a uniform magnetic field B^, the normalized proton wave function is 

exp {—iE^t + ip3yy + ip3zz) fp-iy,p-i^{x), where fpsy,psX^) is the 4- 
component spinor solution |TT|. The form of the spinor in a magnetic field (see Ref. | p.l[ ) 
restricts the analytical evaluation of the neutrino emissivity to fields strong enough so as to 
populate only the ground state for electrons and protons, i.e. t] = rj =0. The only positive 



energy spinor for protons in the chiral representation is then pT| , p!3 



fE*,+ps,\ 




—m 

V J 



Iv=0;P3yix), (6) 
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(pI^ + m*^)^/^ is the effective 



where iV^=o = I / ^2E*{e[+^) , and = 
relativistic Hartree energy. The function Iri=o-p-^y{x) is similar in form as in Ref. fill]. The 
nucleon effective and rest masses are respectively, m* and m = rrin = rrip = 939 MeV. In 
presence of the magnetic field, the wave functions for free electrons ip^iX) have the same 
form as those for protons, but with m* and for protons replaced by the bare mass mg and 
kinetic energy for electrons, respectively. The neutrons and neutrinos/antineutrinos being 
unaffected by Bm, have plane wave functions. 

Using these wave functions, it is straightforward to calculate the transition rate per unit 
volume and is given by 



W 



n2 



EIE^EIE^ V^LyL, 



exp 



{Plx -P2x)'^ + {PSy +P4yf 



X 



2qB^ 

{gv + gAf{pi ■ P2){P3 ■ Pi) + {gv - gAf{pi ■ pa){p3 ■ P2) - {gl - g\)m*^{pi ■ P2) 



X 6{Ei - E2 - E3 - Ei)5{piy - P2y - Psy - PAy)5{Plz " P2z - PZz " PAz 



(7) 



Substituting Eq. (7) in the expression (4) for emissivity, and by the change of variable 
(Psj/ +P4j/) P2,yi the integration over dp^y can be performed to yield a factor qB^nL^. The 
rest of the integrals of Eq. (4) can then be performed in the standard manner . Electron 
capture gives the same emissivity as neutron decay, although in neutrinos, and thus the 
total emissivity (relativistically) for the direct URCA process in nuclear matter (NM) in a 
magnetic field Bm is e^uRCAl-^m) = 25^,(5^) i.e. 



.NM 

'URCA 
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exp 
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PF„ r. 
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pk 



J Pf^Pf, 



(8) 



where fi* = {pf + m*^)^/^ and cos^^u = {p^p^ + p\ — PPp) f^PFnPFe- The threshold factor is 
= 9{pFp + PFe ~ Pf„), where 9{x) = 1 for x > and zero otherwise. For B^ = 0, the 
relativistic expression for the neutrino emissivity from the nucleon direct URCA process is 



.NM 

'URCA 
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(9) 



It was shown [0 that quark matter (QM), if present, the /3-decay (i.e. direct URCA 
process) of d quarks is kinematically allowed through reaction (2) if finite mass (and/or 
quark-quark interaction) is incorporated. The relativistic expression of the neutrino emis- 
sivity for the direct URCA process involving u and d quarks for Bm = and without 
quark-quark interaction is given by R] 



-QM 

'URCA 



{Bm = 0) 



4577r 
840 



Gp COS^ 6c{l - cos 634) fJ'dfJ'ufJ'eT^ 



(10) 



in the usual notation [0]. The emissivity for the (3 decay of s quark for Bm = is similar 
to Eq. (10) with cos^^c replaced by sin6'c. The emissivity for the (3 decay of free d quark 
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in a magnetic field Bm may be obtained from Eq. (8) by substituting gv = Qa = ^ witfi 
fJ'n ~^ f^d, lA, ~^ f^u, and multiplying a color factor 3 for d quark: 



^URCAl^m) = ^^'^i^cos e, [qB^) II - — 1 exp 



2qB„ 



PFuPF, 



Similar expression is obtained for s quark, but is Cabibbo suppressed. The decay of d and 
s quarks is feasible if they satisfy the respective inequality conditions PFu ~ Pf^ ^ Pf^ < 
Pf„ + Pf, and pp^ - Pf, < Pf, < Pf^ + Pf, ■ 

To estimate numerically the various neutrino emissivities for the direct URCA processes 
with and without magnetic field in a neutron star, we describe the nuclear matter and 
electrons within the relativistic Hartree approach in the linear a-uo-p model 0,^. The 
values for the dimensionless coupling constants for the a, uj and p mesons are adopted from 
Ref. which are determined by reproducing the nuclear matter properties at a saturation 



density of uq = 0.16 fm . The variation of magnetic field with density rib from surface to 



center of the star is parametrized by the form [|T^] 

5™K/no) = B^^' + 5o [1 - exp {-(3{n^/n^y]] , (12) 

where the parameters are chosen to be /5 = 10~^ and 7 = 6. The maximum field prevailing 
at the center is taken as i?o = 5 x 10^*^ G and the surface field is B^^^ ~ 10^ G. The number 
of Landau levels populated for a given species is determined by the i?^ and ||11|| . 



In Fig. 1, we show the neutrino emissivity as a function of baryon density at B^ = 
(see Eq. (9)) for the direct URCA process in nuclear matter (denoted by NM) at an interior 
temperature T = 10^ K. Due to momentum conservation, the threshold density at which 
this process occurs, is at rit = 0.346 fm~'^. The variation of emissivity with Ub in presence of 
magnetic field B^ (see Eq. (8)) as seen in the figure may be explained as follows: At very 
low densities Uf, ~ 0.35 — 0.73 the field B^ (as given in Eq. (12)) is rather small ~ 10^^ 
G, and consequently a large number of Landau levels are populated. This gives essentially 
field-free results. At densities Uf, > 0.75 fm~^, the field is strong enough to populate only the 



ground levels of both electrons and protons and would have pronounced quantization 



effects; the critical field for electron is -B.^^*-^-* = 4.414 x 10^'^ G. The emissivity then rapidly 
increases with density and could have values as high as ~ 2 orders of magnitude larger 
than Bm = case at n^, ^ 1.2 fm^^. Hereafter, B^ saturates to a maximum of 5 x 10^*^ G 
so that for > 1.2 fm~'^, higher level states start to populate, and, as in the low density 
situation, results in field-free emissivity values. The central densities Uc of neutron stars with 
maximum masses are also shown in Fig. 1 with (open circles) and without (solid circles) 
the magnetic field. For Bm 7^ star, ric = 1.448 fm~'^ and thus falls above the kernel of 
enhanced emissivity leading to faster cooling compared to the field-free case. 

The neutrino energy losses from direct URCA processes of quark matter composed of 
free u, d and s quarks and e are estimated in the bag model. The current masses of the 
quarks are taken as = 5 MeV, rrid = 10 MeV and rris = 150 MeV, and the bag constant 
as i? = 250 MeV fm~^. In Fig. 1, we display the neutrino emissivity from the /3-decay 
of d and s quarks at B^ = (see Eq. (10)) in quark matter (denoted by QM). The d 
quark /5-decay reactions are kinematically allowed if ~ uq. At densities rif, ~ 0.85 fm~^ 
and above when s quark decay is allowed, the emissivity is increased to about an order of 
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magnitude. This is caused by the large s quark mass which allows the momenta of the free 
particles to deviate appreciably from collinearity which tends to increase the matrix element. 
It was, however, shown that by the inclusion of quark-quark interaction, the neutrino 
emissivities from d and s quark /5-decay are comparable in magnitude. The emissivities 
for the quark direct URCA processes in presence of the magnetic field (see Eq. (11)), 
remain virtually unaltered from the field-free case due to the population of a large number 
of levels in all the quark species. In either case, it is found that £^urca/^urca ~ 10~^. 




FIG. 1. The neutrino emissivities as a function of baryon density from the direct URCA process 
for a magnetic field Bm = (solid line) and for Bm = 10^ — 5 x 10^® G (dashed line) for: nucleons 
in nuclear matter (NM); quarks in quark matter (QM); a nucleon-quark phase transition (NQP); 
nucleons in nuclear matter with hyperons (NM(HY)); a nuclear matter with hyperons to quark 
phase transition (NQP(HY)). The maximum masses of the stars Mmax with these various compo- 
sitions are given for B^ = and those in the parentheses are for Bm = 10^ — 5 x 10^^ G. The 
corresponding central densities are indicated by solid and open circles, respectively. 



In a realistic situation, if quarks at all exist, a star with increasing density from the surface 
to the center would have a pure nucleon phase at the inner crust and core with a possible pure 
quark phase at the center and a mixed nucleon-quark phase (NQP) in between. The mixed 
phase of nucleons and quarks is described following Glendenning The conditions of 

global charge neutrality and baryon number conservation are imposed through the relations 
xQ"' + (1 x)Q'^ = and rib = X^b + (1 ~ x)'^6) where x represents the fractional volume 
occupied by the hadron phase. Furthermore, the mixed phase satisfies the Gibbs' phase 
rules: fip = 2/iu -|- fi^ and P" = P'^. The neutrino energy loss rate in this phase is given 
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t'y ^URCA = X^uRCA + (1 ^ x)^uRCA- The neutrino emissivities for the nucleon-quark phase 
transition are shown in Fig. 1 (denoted by NQP). For Bm = case, with the appearance of 
the quarks at rib = 0.533 fm~^, the emissivity decreases from the corresponding NM case. 
Apart from the reduced emissivity of the quark phase (which being, however small at large 
x), the reduction in the chemical potentials of the nucleons and electrons resulting from 
the requirement of the global charge neutrality and baryon number conservation conditions 
in the mixed phase, primarily causes the decrease in emissivity. For stars with Bm 0, 
the emissivities in the mixed phase are enhanced, particularly in the regime dominated by 
nucleons (i.e. for x > 0.5). The central densities of maximum mass stars fall within the 
mixed phase, and consequently such stars would have faster cooling than pure quark stars. 
The maximum mass NQP stars with and without magnetic field, however, have much smaller 
emissivities than that of the corresponding NM stars, while NQP stars with Bm have nearly 
identical cooling as that of field-free NM stars even though their maximum masses are very 
distinct. 

Since quark matter furnishes both baryon number and negative charge, intuitively, the 
trends exhibited by the emissivities for NQP stars may be anticipated by invoking strange 
baryons, namely hyperons (A's, S's and S's). The /3-equilibrium conditions then generalize 
to fii = hiHn — Qifie, where bi and qi are the baryon number and charge for the ith particle. 
Since the hyperons are more massive than the protons, the effect of the magnetic field on 
their direct URCA processes is negligible. Because of the large uncertainties in the hyperon- 
nucleon interactions even at nuclear density, for a conservative estimate of the emissivities 
we set the nucleon-meson and hyperon-meson coupling constants equal. Furthermore, the 
critical density for nucleon direct URCA process is nearly identical to the hyperon threshold 
density in the relativistic mean field model, and the emissivities from the hyperon direct 
URCA processes are about 5-100 times less than that from the nucleons Therefore, we 
shall present emissivity vs rify results only for the nucleon direct URCA process in presence 
of hyperons. This is shown in Fig. 1 and denoted by NM(HY). With the appearance of 
hyperons, the reduction in the chemical potentials of the nucleons and electrons required 
by the baryon number conservation and charge neutrality condition causes a substantial 
reduction of the emissivity compared to that from NM. In fact, with increasing density 
when hyperon abundances grow rapidly, the emissivities gradually decrease. For Bm 7^ 0, 
only the ground Landau levels for e and p are populated over a considerable density range 
in this matter. Consequently, the emissivities with Bm in NM(HY) stars are significantly 
larger than that for the corresponding field-free stars. 

Allowing now baryon to quark phase transition, the emissivity displayed in Fig. 1 (de- 
noted by NQP(HY)) for Bm = is larger than that from the NQP matter. This is caused 
by the delayed appearance of quarks in hyperon rich matter, so that the the total emissivity 
is primarily dominated by the nucleons. In presence of the field, the total emissivity of 
NQP(HY) matter is about an order of magnitude larger for the maximum mass star and 
therefore leads to faster cooling compared to the corresponding field-free star. 

Within the non-relativistic (but interacting) approximation for the specific heat Cy and 
neutrino emissivity, the time for the center of a NM star to cool by the direct URCA process 
to a temperature Tg at Bm = may be estimated to be At = — J{cy/e^^Qji^)dT ~ lOTg"^ s. 
In contrast, for Bm 7^ 0, the NM star's center cools faster with At ~ 0.5Tg~^ s. By invoking 
quarks and/or hyperons, the decrease in emissivity is much more compared to that of the 
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specific heat resulting in slow cooling of the stars center. The typical time scale associated 
with the propagation of thermal signals through the outer core and the crust to the surface 
before the sudden temperature drop is quite high ~ 1 to 100 yr, depending on the crustal 
composition and relative sizes of the crust and the core and thus upon the equation of state. 
Therefore, it seems to be quite difficult to distinguish observationally from the effects of 
direct URCA process, the interior constitution of a star. 

Throughout our discussion we have assumed that the electron is the only lepton. If the 
triangle inequalitj pp^+pp^ > pp^ is satisfied then nucleon direct URCA process with muons 
will occur; the threshold density for this process is higher than electrons since > nif.. The 
/3-equilibrium condition = /^^ moreover implies that the emissivity for URCA process with 
muons is same as that for the corresponding process with electrons. In the present model the 
nucleon direct URCA processes are not permitted at densities rih < 0.34 fm~^. In this outer 
core region, the dominant neutrino emission process are then the modified URCA processes 
of reaction (1) for which the emissivity is smaller by a factor ~ (T/Tp)^ than the direct 
URCA processes. 

At certain densities and temperature T < Tc ^ 10^ — lO^'' K, the nucleon superfiu- 
idity may set in. The specific heat and direct URCA rate are then reduced by a factor 
~ exp(— A/T), where A is the larger of the neutron and proton gaps. The modified URCA 
rates are, however, reduced by a factor ~ exp(— 2A/T). In presence of a magnetic field, the 
superfluid protons are believed to form a type II superconductor in the outer core within the 
density range 0.7no < rib < 2no, and the estimated lower and upper critical magnetic fields 
are respectively Hd ~ 10^^ G and Hc2 ~ 3 x 10^^ G fl^. With the choice of variation of 



Bm with rib (see Eq. (12)), the field is 10^^ — 10^^ G at the bulk of the outer core and could 
form a superconducting region at T < T^, while the inner core and center with Bm ~ 10^^ 
G is in the normal state without superconductivity. 

In conclusion, the neutrino emissivities for all the cases studied here are found to be 
dramatically enhanced in a magnetic field compared to that from the non-magnetized stars. 
However, for certain stars unambiguous determination of the interior constituents may be 
difficult. There can be stars with the same composition, as for example the NM stars in a 
magnetic field but with slightly different masses of 1.55Mq and 1.60Mq having their central 
densities at 0.72 fm~^ and 0.88 fm~^ residing below and within the kernel of fast cooling 
respectively, and thus have completely different emissivities. On the other hand, a NM 
star with Bm = and a NQP star in a magnetic field, though possessing different interior 
compositions, have nearly identical emissivities. It is also found that when pure nuclear 
matter is injected with nonleptonic negative charges, namely hyperons and quarks, the 
emissivities turn out to be smaller than that from the nuclear matter. It has been already 
demonstrated [|1^ that nonleptonic negative charges cause a softening of the equation of 
state. We thus arrive at a general result that when matter contains nonleptonic negative 
charges, the maximum masses of the stars are smaller with a suppression of the neutrino 
emissivity than that of the pure nuclear matter with and without a magnetic field. 
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